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PREFACE 


|  In  August  1983,  the  US  Army  Engineer  Waterways  Experiment  Station  (WES) 

was  requested  by  the  Norfolk  District  Corps  of  Engineers  to  conduct  an  in¬ 
vestigation  of  the  possible  sedimentation  changes  in  the  federally  maintained 
channels,  general  sedimentation  in  the  lower  James  River,  and  changes  in 
|  flushing  near  Craney  Island  caused  by  the  proposed  1-664  Bridge-Tunnel  com¬ 

plex.  This  study  was  funded  by  the  Federal  Highway  Administration  and  the 
Virginia  Department  of  Highways  and  Transportation. 

The  study  was  conducted  from  August  1983  to  August  1984  by  personnel  of 
|  the  Hydraulics  Laboratory,  WES,  under  the  general  direction  of  Messrs.  H.  B. 

Simmons  and  F.  A.  Herrmann,  Jr.,  former  and  present  Chiefs  of  the  Hydraulics 
Laboratory,  respectively;  R.  A.  Sager,  Chief  of  the  Estuaries  Division;  E.  C. 
McNair,  Chief  of  the  Sedimentation  Branch;  and  R.  A.  Boland,  Chief  of  the  Hy¬ 
drodynamics  Branch.  The  project  was  conducted  by  Messrs.  S.  B.  Heltzel, 

M.  J.  Trawle,  and  R.  F.  Athow,  Estuaries  Division.  Mr.  Heltzel  prepared  this 
report.  Mr.  D.  Stewart,  Estuaries  Division,  also  participated  in  this  study. 
This  report  was  edited  by  Information  Technology  Laboratory  personnel 
Mrs.  Beth  Burris  and  Mrs.  Marsha  Gay  with  Mrs.  Chris  Habeeb  coordinating  the 
final  layout. 

Mr.  J.  R.  Melchor,  Norfolk  District  Corps  of  Engineers,  made  valuable 
contributions  to  this  project. 

Col  Dwayne  G,  Lee,  EN,  is  the  Commander  and  Director  of  WES. 

Dr.  Robert  W.  Whalin  is  the  Technical  Director. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Multiply 

By 

To  Obtain 

acres 

0.4047 

hectares 

cubic  feet 

0.02831685 

cubic  metres 

degrees  (angle) 

0.01745329 

radians 

feet 

0.3048 

metres 

miles  (US  nautical) 

1.852 

kilometres 

miles  (US  statute) 

1.609344 

kilometres 

square  miles  (US  statute) 

2.589998 

square  kilometres 
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1-664  BRIDGE-TUNNEL  STUDY,  VIRGINIA 
SEDIMENTATION  AND  CIRCULATION  INVESTIGATION 


PART  I :  INTRODUCTION 

The  James  River  Estuary 


Background 

1.  The  James  River  is  a  narrow,  funnel-shaped,  shallow  estuary  with  a 
mean  low-water  average  depth  of  12. I  ft.*  Riverbed  channels  have  an  average 
depth  of  32.8  ft.  Some  channel  scour  holes  exceed  98.4  ft  in  depth.  At 
Hampton  Roads,  where  the  James  River  becomes  wider  as  the  Elizabeth  and 
Nansemond  Rivers  enter  the  estuary,  the  harbor  deepens  to  about  59  ft  with  a 
channel  dredged  at  its  mouth.  Figure  1  details  the  bathymetry  of  the  James 
River  estuary. 

2.  The  James  River,  a  partially  mixed  estuary  of  the  Chesapeake  Bay 
estuary  system,  drains  the  geological  provinces  of  the  Blue  Ridge,  the  Pied¬ 
mont,  and  the  Coastal  Plain  of  Virginia.  The  three  major  tributaries  entering 
the  James  River  are  the  north  branch  of  the  James  River,  with  a  drainage  area 
of  6,757  square  miles;  the  Appomattox  River,  which  enters  the  James  River  at 
Hopewell,  VA,  with  a  drainage  area  of  1,344  square  miles;  and  the  Chickahominy 
River,  which  enters  the  James  River  along  Hog  Island,  with  a  drainage  area  of 
248  square  miles  (Nichols  1972).  Minor  tributaries  Include  the  Pagan,  the 
Nansemond,  the  Elizabeth,  and  the  Warwick  Rivers.  The  three  major  tributaries 
have  gaging  stations  and  have  recorded  the  following  average  discharges:  the 
North  James  River,  7,098  cfs;  the  Appomattox  River,  1,095  cfs;  and  the 
Chickahominy  River,  247  cfs  (Nichols  1972).  Minor  tributaries  have  not  been 
gaged  and  are  assumed  to  have  comparatively  low  discharges. 

Circulation** 

3.  During  a  mean  tide,  tidal  currents  vary  from  -1.15  to  +1.48  fps  from 
Newport  News  to  Jamestown.  Tidal  currents  reach  their  peak  3  hr  after  slack. 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  found  on  page  3. 

**  Summarized  from  Nichols  (1972), 


Figure  1.  Bathymetric  chart  of  the  James  River  estuary,  based  on  US  Coast 
and  Geodetic  Survey  charts  400  and  529.  BP  is  Brandon  Point;  RL,  Rockland¬ 
ing  Shoal  Channel;  WR,  Wreck  Shoal;  BB,  Burwell  Bay;  MP,  Mulberry  Point; 
OPC,  Old  Point  Comfort;  RAP,  Rappahannock  River;  e,  ebb  channel;  f,  flood 
channel  (Nichols  1972)  (Originally  published  by  the  Geological  Society  of 
America  as  "Sediments  of  the  James  River  Estuary,  Virginia,"  by  Maynard  M. 
Nichols,  in  Environmental  Framework  of  Coastal  Plain  Estuaries,  Geological 
Society  of  America  Memoir  133,  pp  169-212) 


There  is  a  density  circulation  in  the  James  River  that  is  potentially  impor¬ 
tant  in  transporting  sediment  up  the  estuary. 

4.  Salinity  varies  with  distance  downstream  from  0  ppt  at  Richmond,  VA, 

r0  to  24  ppt  at  the  mouth.  High  river  inflow  produces  the  greatest  stratifica- 

tion  in  the  upper  estuary.  The  limit  of  saltwater  intrusion  varies  from  Hog 
Island,  24  miles  above  the  mouth,  during  high  river  discharge  to  a  location 
approximately  54  miles  above  the  mouth  during  low  discharge  in  the  summer  and 
:•  fall. 

5.  The  seasonal  variations  in  river  inflow  change  the  estuary  from 
moderately  stratified  to  well  mixed.  At  low  river  inflow,  tidal  mixing  re¬ 
duces  stratification,  producing  a  vertically  homogeneous  estuary.  Figure  2 
shows  the  distribution  of  surface  salinity  averaged  over  one  to  two  tidal 
cycles  or  at  slack  water. 

5 


SURFACE  SALINITY,  %. 


80  70  60  -• — Dl  STANCE  UPSTREAM  20  10  0  km 


TOTAL  SUSPENDED  SEDIMENTS,  mg/liUr 

Figure  2.  Distribution  of  total  suspended  sediment  con¬ 
centrations  and  surface  salinity  along  the  estuary  chan¬ 
nel,  based  on  average  values  over  one  to  two  tidal  cycles 
between  31  miles  and  the  mouth,  and  on  slack-water  con¬ 
centrations  between  31  and  50  miles  above  mouth, 

11-20  March  1965  (Nichols  1972)  (Originally  published  by 
the  Geological  Society  of  America  as  "Sediments  of  the 
James  River  Estuary,  Virginia,"  by  Maynard  M.  Nichols,  in 
Environmental  Framework  of  Coastal  Plain  Estuaries,  Geo¬ 
logical  Society  of  America  Memoir  133,  pp  169-212) 


Sediment  character- 
istics  and  distribution 

6.  Suspended  sediment  concentrations  increase  from  an  average  of  about 
30  mg/ £  above  Jamestown  to  more  than  100  mg/ £  at  the  limit  of  saltwater  intru¬ 
sion  (Deep  Water  Shoal).  Concentrations  decrease  to  12  mg/£  at  the  mouth. 
Although  concentrations  vary  during  the  tidal  cycle,  this  overall  pattern  per¬ 
sists.  Measured  concentrations  are  low  compared  with  other  east  coast  estu¬ 
aries.  Figure  2  shows  the  distribution  of  total  suspended  sediment  concentra¬ 
tion  averaged  over  one  to  two  tidal  cycles  or  at  slack  water. 

7.  The  turbidity  maximum  is  located  at  the  upstream  limit  of  salinity 
intrusion.  This  maximum  persists  most  of  the  year,  being  most  pronounced  in 
the  spring  when  river  inflow  is  high  and  weak  in  the  fall  when  inflow  is  low. 
Nichols  (1972)  reports  the  following:  "Located  close  to  or  slightly  upstream 
of  the  0.5  ppt  isohaline,  the  maximum  shifts  upstream  with  landward  penetra¬ 
tion  of  salt  water  from  spring  to  fall.  This  trend  indicates  that  the  posi¬ 
tion  and  intensity  of  the  maximum  may  be  partly  controlled  by  the  inner  limit 
of  salt  water  and  in  turn  by  the  magnitude  of  river  inflow." 


8.  Nichols  (1972)  also  states:  "It  is  evident  that  about  20  percent  of 
the  suspended  load  remains  in  suspension  at  slack  water,  whereas  the  greater 
part  is  alternately  suspended  and  settled  out.  These  changes  are  not  sym¬ 
metrical  about  slack  water  or  maximum  current.  Concentrations  near  the  sur¬ 
face  are  slightly  higher  at  the  end  of  ebb  than  at  the  end  of  flood,  a  trend 
ascribed  to  shifting  of  the  longitudinal  concentration  gradient."  Figure  3 
shows  patterns  of  deposition  and  erosion  in  the  James  River  estuary. 
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Figure  3.  Patterns  of  deposition  and  erosion  in  James  River  estuary  deter¬ 
mined  from  depth  changes  over  a  70-year  period.  Prominent  zone  of  deposi¬ 
tion  (a)  in  the  lower  estuary;  Burwell  Bay  (b) ,  the  main  site  of  deposition. 
Dotted  line  is  a  13-ft  depth  contour.  Arrows  represent  direction  of  net 
flow  along  bottom  (Nichols  1972)  (Originally  published  by  the  Geological 
Society  of  America  as  "Sediments  of  the  James  River  Estuary,  Virginia,"  by 
Maynard  M.  Nichols,  in  Environmental  Framework  of  Coastal  Plain  Estuaries, 
Geological  Society  of  America  Memoir  133,  pp  169-212) 


Scope 


9.  The  navigation  channel  sedimentation  was  to  be  studied  using  the 
Elizabeth  River/lower  James  River  numerical  mesh  developed  for  the  Norfolk 
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Harbor  and  channels  deepening  study  (Berger  et  al.  1985).  A  new  mesh  was 
developed  for  the  general  sedimentation  study.  Data  for  the  circulation  study 
were  provided  by  the  Virginia  Institute  of  Marine  Science  (VIMS)  (Fang  et  al. 
1972)  from  a  1972  VIMS  study  of  the  proposed  bridge-tunnel  crossing  conducted 
with  the  James  River  physical  model  at  the  US  Army  Engineer  Waterways 
Experiment  Station  (WES). 


Proposed  Bridge-Tunnel  Complex 

10.  The  proposed  1-664  Bridge-Tunnel  complex  consists  of  (a)  two  man¬ 
made  islands,  (b)  a  tunnel  section  connecting  the  two  islands,  and  (c)  a 
bridge  portion  connecting  the  south  island  and  the  south  shore.  Figure  4 
shows  the  location  of  the  bridge-tunnel  complex. 

Purpose 

11.  The  purpose  of  this  study  was  to  determine  the  effects  of  the  pro¬ 
posed  1-664  James  River  Bridge-Tunnel  complex  on  (a)  sedimentation  in  the 
federally  maintained  channels  (Newport  News,  Norfolk  Harbor,  and  Elizabeth 
River) ,  (b)  general  sedimentation  in  the  lower  James  River  (specifically  those 
regions  containing  oyster  and  clamming  beds),  (c)  changes  in  overall  flushing 
characteristics,  and  (d)  changes  in  current  velocities  and  flushing  near  the 
Craney  Island  disposal  site.  Items  (a)  and  (b)  were  to  be  numerical  model 
studies  and  items  (c)  and  (d)  were  to  be  desk  studies. 

12.  Tests  involving  the  impact  of  the  proposed  1-664  Bridge-Tunnel  com¬ 
plex  on  wave  conditions  in  the  Newport  News  Small-Boat  Harbor  were  investi¬ 
gated  at  WES  and  are  reported  in  Bottin  (1984). 
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Figure  4.  1-664  crossing  of  Hampton  Roads 
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PART  II:  THE  MODELS 


James  River  Physical  Model 


13.  The  1972  VIMS  study  (Fang  et  al.  1972)  was  conducted  in  the  exist¬ 
ing  James  River  model,  which  reproduces  the  entire  James  River  estuary,  a  por¬ 
tion  of  lower  Chesapeake  Bay,  and  about  200  square  miles  of  the  Atlantic  Ocean 
as  shown  in  Figure  5.  The  tidal  portions  of  all  major  tributary  streams  were 
reproduced,  including  the  Elizabeth,  Nansemond,  Pagan,  Warwick,  Chickahominy , 
and  Appomattox  Rivers.  The  model  is  about  550  ft  long  and  130  ft  wide  at  the 
widest  point. 

14.  Since  gravitational  forces  are  predominant  in  tidal  flows,  it  can 
be  determined  that  the  model  and  prototype  Froude  numbers  F  must  be  equal. 
Therefore 


(1) 


where 

2 

V  /(gL)  *  Froude  number 

V  =  velocity  of  flow 
g  «*  acceleration  of  gravity 
L  =  length  dimension 

m,p  -  subscripts  indicating  model  and  prototype,  respectively 
Geometric  scales  of  the  model  are  1:1,000  horizontally  and  1:100  vertically, 
reflecting  a  geometric  distortion  ratio  of  10:1.  These  dimensions  and 
Froude  model  laws  defined  the  following  model-to-prototype  ratios: 


Characteristic 

Time 

Velocity 

Volume 

Discharge 


Scale  Relations 
Model  Prototype 

1:100 

1:10 

1:100,000,000 

1:1,000,000 


The  model-to-prototype  ratio  for  salinity  is  1:1. 


Figure  5.  James  River  physical  model  limits 


15.  The  physical  model  of  Chesapeake  Bay  is  located  on  Kent  Island  in 
Matapeake,  MD.  It  is  an  8.6-acre,  fixed-bed  model.  The  area  reproduced  in 
the  model  extends  from  approximately  30  miles  offshore  in  the  Atlantic  Ocean 
to  the  heads  of  tide  for  all  tributaries  emptying  into  the  Chesapeake.  The 
entire  length  of  the  Chesapeake  and  Delaware  (C&D)  Canal  and  a  portion  of 
Delaware  Bay  are  also  modeled.  Overbank  geometry  is  reproduced  to  the  +20-ft 
contour.  Model  limits  are  shown  in  Figure  6. 

16.  The  model-to-prototype  scale  ratios  are  the  same  for  this  model  as 
for  the  model  of  the  James  River. 

The  Numerical  Models 

17.  The  finite  element  hydrodynamic  model  RMA-2V  was  used  in  both  the 
navigation  channel  sedimentation  study  and  the  general  sedimentation  study. 

The  finite  element  transport  model  STUDH,  operating  in  a  noncohesive  sediment 
mode,  was  used  in  the  navigation  channel  sedimentation  study  to  evaluate 
changes  in  shoaling  in  the  federally  maintained  channels.  The  cohesive  sedi¬ 
ment  mode  was  used  in  the  navigation  channel  sedimentation  study  to  evaluate 
changes  in  shoaling  in  the  Elizabeth  River  and  in  the  general  sedimentation 
study.  RMA-2V  and  STUDH  are  components  of  the  TABS-2  numerical  modeling  sys¬ 
tem,  the  following  description  of  which  was  taken  from  McAnally  et  al.  (1983) 
and  Thomas  and  McAnally  (1985). 

Finite  element  mesh  generation 

18.  The  finite  element  models  RMA-2V  and  STUDH  require  that  a  repre¬ 
sentation  of  the  area  to  be  studied  be  put  in  a  digital  form.  This  represen¬ 
tation  is  in  the  form  of  a  set  of  computational  points  that  have  been  located 
with  respect  to  a  reference  system.  These  computational  points,  called  nodes, 
are  assigned  sequential  numbers  and  bed  elevations.  The  nodes  are  connected 
to  one  another  by  lines  that  create  either  triangular  or  quadrilateral  ele¬ 
ments.  The  nodes  midway  between  corners  are  called  midside  nodes. 

The  hydrodynamic  model,  RMA-2V 

19.  The  hydrodynamic  model,  RMA-2V,  solves  the  depth-integrated  equa¬ 
tions  of  fluid  mass  and  momentum  conservation  in  two  horizontal  directions. 

The  present  model,  an  improvement  of  an  earlier  version,  RMA-2  (Norton  and 


King  1977) ,  is  formulated  in  terms  of  velocities  and  turbulent  exchange 
coefficients. 
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20.  The  finite  element  method  using  the  method  of  weighted  residuals  is 
used  to  solve  the  conservation  of  mass  and  momentum  equations.  Individual 
elements  may  be  either  quadrilaterals  or  triangles  and  may  have  curved  (para¬ 
bolic)  sides.  Shape  functions  are  quadratic  for  flow  and  linear  for  depth. 
Integration  in  space  is  performed  by  Gaussian  integration.  Derivatives  in 
time  are  replaced  by  a  nonlinear  finite  difference  approximation. 

21.  The  finite  element  solution  is  fully  implicit  and  a  set  of  simul¬ 
taneous  equations  is  solved  by  Newton-Raphson  iteration.  The  solution  is 
achieved  using  a  front-type  matrix  solver  that  assembles  a  portion  of  the 
matrix  and  solves  that  portion  before  assembling  the  next  portion  of  the  ma¬ 
trix.  The  front  solver's  efficiency  is  largely  independent  of  bandwidth  and 
thus  does  not  require  as  much  care  in  formation  of  the  computational  mesh  as 
do  traditional  solvers. 

The  sediment  transport  model,  STUDH 

22.  Convection-diffusion  equation.  The  sediment  transport  model, 

STUDH,  solves  the  depth-integrated  convection-diffusion  equation  in  two  hori¬ 
zontal  dimensions  for  a  single  sediment  constituent.  The  form  of  the  solved 
equation  is 


2£  + 
at 


u 


|£  + 

3x 


ac 

3z 


+  c^C  +  <*2 


(2) 


where 

3 

C  -  concentration  of  sediment,  kg/m 
t  *  time,  sec 

u  *  flow  velocity  in  x-direction,  m/sec 
x  *  primary  flow  direction,  m 
w  *  flow  velocity  in  z-direction,  m/sec 
z  ■  direction  perpendicular  to  x  ,  m 

2 

D  *  effective  diffusion  coefficient  in  x-direction,  m  /sec 

X  2 
D  »  effective  diffusion  coefficient  in  z-direction,  m  /sec 
z 

a.  *  coefficient  for  concentration-dependent  source/sink  term,  1/sec 

1  3 

02  ■  coefficient  of  source/sink  term,  kg/m  /sec 

23.  Bed  shear  stress.  The  bed  shear  stress,  r,  ,  takes  the  form: 
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p  *  water  density 
u^  =  shear  velocity 

24.  Manning  shear  stress  equation.  The  Manning  form  of  the  shear 
stress  equation  was  used  in  this  study: 


CME  f  T>‘ 


where 

u  *  flow  velocity 

n  -  Manning' 8  roughness  coefficient 
g  *»  acceleration  due  to  gravity 
CME  =  coefficient  of  1  for  SI  units  and  1.486  for  non-Sl  units 
D  *  flow  depth 

25.  Sand  transport.  This  supply  of  sediment  to  and  from  the  bed  for 
noncohesive  bed  material  (sand)  is  controlled  by  the  transport  potential  of 
the  flow  and  the  availability  of  material  in  the  bed.  The  bed  source  term  is 
calculated  by 


Ceq  -  C 
t 


where 


S  *  source  term 


Ceq  -  equilibrium  concentration 
tc  ■  characteristic  time  for  effecting  the  transition 

26.  The  value  of  Ceq  can  be  determined  from  any  of  several  transport 
relations.  The  sand  version  of  STUDH  uses  the  Ackers-White  formula  (1973), 
which  performed  satisfactorily  in  tests  by  WES  and  others  (White,  Milll,  and 
Crabbe  1975;  Swart  1976). 

27.  Clay  transport.  Deposition  rates  for  clay  beds  were  calculated 
with  the  equations  of  Krone  (1962): 


C  [1  -  —  )  for  C  <  C 

°  Tdj 


P  c5l3(i-Ti)  for  c>cc 


where 

Vg  »  fall  velocity  of  a  single  particle 

t,  =  critical  shear  stress  for  deposition 
d 

C  *  critical  concentration  *  300  mg/ l 

V  =  V  /C  4/3 
k  sc 

28.  Erosion  rates  were  computed  by  Ariathurai's  equation  (Ariathurai, 
MacArthur,  and  Krone  1977)  for  particle-by-particle  erosion.  The  source  term 
is  computed  by 


-Hi-1 


where 


P  *  erosion  rate  constant 

«=  critical  shear  stress  for  particle  erosion 


29.  STUDH  is  a  descendent  of  the  model  SEDIMENT  II  (Ariathurai, 
MacArthur,  and  Krone  1977)  developed  under  the  direction  of  R.  B.  Krone  at  the 
University  of  California,  Davis. 
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PART  III:  MODELING  PROCEDURES 


Circulation 

30.  The  primary  source  of  information  for  this  portion  of  the  study  is 
the  1972  VIMS  study  (Fang  et  al.  1972)  on  the  proposed  bridge-tunnel  crossing. 
The  study  was  conducted  in  the  James  River  hydraulic  model  at  WES  from  15  to 
30  May  1972.  The  purpose  of  the  hydraulic  model  study  was  to  determine  the 
effects  of  the  proposed  1-664  river  crossing  structures  on  the  tides,  cur¬ 
rents,  and  distribution  of  salinity  and  sediments  in  the  reach  between  Old 
Point  Comfort  and  the  existing  James  River  Bridge.  In  addition  to  the  exist¬ 
ing  (base)  conditions,  two  basic  configurations  were  considered.  The  first 
was  the  base  configuration  plus  the  new  islands  for  the  1-664  Hampton  Roads 
Bridge-Tunnel  crossing  (Plan  1A)  and  the  second  James  River  Bridge  crossing. 
The  second  configuration  was  identical  with  the  first,  but  added  a  westward 
triangular  extension  of  the  Craney  Island  disposal  area  (Plan  IB) . 

Physical  model  operating  conditions 

31.  Throughout  the  series  of  model  tests,  the  freshwater  inflow  at 
Richmond,  VA,  was  maintained  at  7,500  cfs,  which  is  the  average  yearly  fresh¬ 
water  flow.  The  total  freshwater  inflow  from  the  Appomattox,  Chickahominy, 
and  other  tributaries  was  2,000  cfs.  The  model  was  operated  with  a  mean  tide, 
and  the  tidal  gage  near  Thimble  Shoal  was  used  as  the  control  gage.  The 
Chesapeake  Bay  salinity  was  maintained  at  24.2  ppt.  These  particular  boundary 
conditions  correspond  to  a  Condition  4  run  in  the  Chesapeake  Bay  physical 
model  for  the  Norfolk  Harbor  and  channels  deepening  study  (Richards  and  Morton 
1983)  and  can  be  compared  to  those  data. 

Data  collected 

32.  Table  1  summarizes  what  data  were  obtained  by  VIMS  for  the  base  and 
Plan  1A.  These  data  included  velocity  and  salinity  measurements  taken  at  sur¬ 
face  (S),  middepth  (M) ,  and  bottom  (B).  Figures  7-10  show  the  locations  of 
the  stations  sampled  with  the  number  of  depths  sampled.  Locations  of  the 
tidal  height  stations  are  also  shown.  All  data  were  collected  half-hourly 
over  a  tidal  cycle. 

Data  collection  procedure 

33.  Tidal  heights  were  measured  for  three  tidal  cycles  at  three  tide 
gage  stations:  Hampton  Roads,  which  is  located  in  Norfolk  at  the  Navy 


Current  measuring  station  Figure  8.  Current  measuring  station 

locations  (base)  locations  (Plan  1A) 


Figure  9.  Salinity  measuring  station  Figure  10.  Salinity  measuring  station 

locations  (base)  locations  (Plan  1A) 


Shipyard  on  the  Elizabeth  River;  Newport  News  Shipbuilding  and  Dry  Dock  in 
Newport  News;  and  Miles,  which  is  located  just  upstream  from  the  James  River 
Bridge  on  the  Newport  News  side  of  the  channel. 

34.  In  the  VIMS  study,  a  simple  technique  was  devised  to  measure  the 
current  directions.  Direction  rosettes  were  placed  on  the  model  floor  under¬ 
neath  the  current  meters  and  a  piece  of  thread  was  attached  to  the  rear  of  the 
meter  bracket.  The  thread  trailing  behind  the  meter  indicated  the  current 
direction.  Fang  et  al.  (1972)  indicate  that  at  most  stations  where  the  cur¬ 
rent  direction  did  not  vary  widely  and  frequently,  the  direction  reading  could 
be  accurate  to  ±15  deg. 

35.  Surface  current  photographs  were  also  taken  consisting  of  time 
exposures  of  confetti  moving  along  the  water  surface.  Actual  current  magni¬ 
tudes  can  be  measured  from  the  photographs,  which  were  taken  hourly  during  the 
tidal  cycle.  Figure  11  shows  the  area  over  which  the  surface  current  photo¬ 
graphs  were  taken. 

Analysis  of  data 

36.  To  analyze  the  surface  current  photographs,  a  sample  grid  (Fig¬ 
ure  12)  was  established  to  allow  for  consistent  reading  of  point  values.  This 
grid  was  then  placed  over  each  photograph,  and  the  current  magnitudes  were 
determined  for  each  location  and  compared  with  the  base  conditions  (before  the 
structures  were  included  in  the  model) . 

37.  At  Range  3  (Figure  7),  hourly  discharge  calculations  were  made  to 
determine  the  changes  in  the  tidal  prism  and  residence  time  upstream  of  the 
range.  Base  and  plan  data  at  other  stations  were  compared  to  the  Range  3  data 
for  an  estimate  of  changes  caused  by  the  bridge. 

Navigation  Channel  Sedimentation 


Grid  evolution 

38.  The  numerical  model  mesh  to  study  the  navigation  channel  sedimenta¬ 
tion  was  a  modified  version  of  the  mesh  used  to  study  Elizabeth  River  shoaling 
for  the  Norfolk  Harbor  and  channels  deepening  study  (Berger  et  al.  1985).  The 
Elizabeth  River  mesh  contained  1,496  computational  nodes  and  407  elements. 

This  mesh  was  revised  to  include  a  representation  of  the  1-664  Bridge-Tunnel 
complex.  Details  for  the  alignment  of  the  bridge/tunnel  used  in  this  mesh 
were  provided  by  the  Virginia  Department  of  Transportation.  This  alignment 
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Figure  11.  Surface  current  photograph  coverage 
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Figure  12.  Surface  current  photograph  comparison  grid,  48  locations 
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1,582  computational  nodes  and  436  elements.  Figure  13  is  the  numerical  model 
mesh  used  in  this  portion  of  the  study. 

Hydrodynamic  modeling  procedures 

39.  The  hydrodynamic  base  condition  runs  developed  for  the  mesh  used  to 
study  Elizabeth  River  shoaling  for  the  Norfolk  Harbor  and  channels  deepening 
study  (Berger  et  al.  1985)  were  used  as  the  base  conditions  to  study  sedimen¬ 
tation  in  the  navigation  channels.  To  simulate  the  1-664  Bridge  and  islands, 
the  friction  coefficient  was  adjusted  for  those  elements  that  represented  the 
islands  and  bridge.  This  change  caused  the  two  islands  to  block  the  flow  of 
water  and  the  bridge  to  restrict  the  flow  of  water.  Numerical  model  parame¬ 
ters  established  for  the  original  mesh  at  other  locations  were  not  changed. 
Therefore  it  was  not  necessary  to  reverify  this  mesh. 

40.  The  hydrodynamic  plan  condition  runs  were  then  made  using  the  modi¬ 
fied  grid.  The  four  sets  of  boundary  conditions  established  for  the  hydrody¬ 
namic  base  condition  runs  (see  paragraphs  51-53)  were  used  for  the  plan  runs. 
Sedimentation  modeling  procedures 

41.  Cohesive.  To  determine  the  effects  of  the  1-664  Bridge-Tunnel 
complex  on  the  shoaling  in  the  Elizabeth  River  and  Norfolk  Harbor  Channels, 
the  sediment  model  was  used  with  cohesive  sediment  calculations.  Since  the 
base  condition  had  been  run  for  the  Norfolk  Harbor  and  channels  deepening 
study  (Berger  et  al.  1985),  those  results  were  used  directly.  For  the  plan 
runs,  the  hydrodynamic  results  for  the  revised  mesh  were  used.  The  various 
model  parameters  established  during  adjustment  of  the  sediment  model  used  for 
the  Norfolk  Harhor  and  channels  deepening  study  were  used  for  these  plan  runs 
This  included  sediment  concentration  boundary  conditions.  Results  from  these 
runs  were  compared  with  the  Norfolk  Harbor  and  channels  deepening  study  plan 
runs  and  the  verification  was  found  to  be  intact. 

42.  Noncoheslve.  To  determine  the  effects  of  the  I-6o4  Bridge-Tunnel 
complex  on  shoaling  in  the  Newport  News  Channel,  the  sediment  model  was  used 
with  noncohesive  sediment  calculations.  The  sediment  model  was  run  with  the 
four  sets  of  base  and  plan  hydrodynamic  input. 

43.  The  length  of  the  sedimentation  model  time-step  selected  was 
1,800  sec  (0.5  hr).  This  was  small  enough  to  prevent  advection  of  sediment 
completely  through  an  element  during  one  time-step.  The  computational  mesh 
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was  identical  with  that  of  the  hydrodynamic  model.  Other  model  parameters 


were  generally  as  follows: 

_ Parameter _  Value 

Crank-Nicholson  implicitness  factor  0.70 

Manning's  n  roughness  0.017 

Effective  particle  diameter  for  transport,  mm  0.20 

Effective  settling  velocity,  m/sec  0.01 

3 

Boundary  concentration,  kg/m  0.008 


44.  This  modeling  effort  was  conducted  assuming  noncohesive  sediment 
transport  with  the  given  effective  particle  diameter  for  transport  as  being 
coarse  silt  to  very  fine  sand. 

45.  Analysis  of  prototype  surveys  during  the  period  1970-1977  indicated 
an  average  shoaling  rate  of  about  0.04  ft/year  in  the  Newport  News  Channel, 
i.e.,  a  nearly  stable  channel.  After  minor  adjustments,  the  noncohesive 
(sand)  run  of  the  numerical  sediment  model,  STUDH,  agreed  well  with  the  ob¬ 
served  shoaling  rate  of  about  0.04  ft/year. 

46.  Results  from  the  four  boundary  condition  runs  were  combined  as 
described  in  Berger  et  al.  (1985)  and  compared  with  the  base.  The  sediment 
results  from  plans  1-4  were  combined  using  the  percentages  in  paragraph  66  and 
compared  to  the  base  results  combined  in  the  same  way. 


General  Sedimentation 

47.  The  lower  James  River  contains  valuable  oyster  beds  that  could 
be  affected  by  the  bridge-tunnel  complex.  Figure  14  shows  the  location  of 
these  oyster  beds  and  clamming  grounds.  This  figure  is  a  composite  of  various 
location  maps  provided  by  the  Norfolk  District.  The  general  sedimentation 
study  was  performed  to  specifically  evaluate  the  sedimentation  changes  within 
the  confines  of  the  oyster  beds  and  clamming  grounds.  Locations  4,  5,  and  8 
were  not  evaluated  since  these  areas  are  not  actively  harvested. 

Development  of  the  numerical  model  mesh 

48.  The  numerical  model  mesh  was  developed  in  two  stages.  The  first 
mesh  incorporated  the  general  features  of  the  study  area.  In  the  final  mesh, 
the  oyster  grounds  were  resolved  in  greater  detail  and  the  bridge-tunnel  was 
more  accurately  represented.  The  Nansemond  River  was  crudely  represented  but 
did  simulate  the  appropriate  tidal  prism  volume.  The  final  mesh  contained 


1,493  nodes  and  464  elements.  Figure  15  shows  the  completed  numerical  mesh 
for  this  study. 

49.  Two  National  Ocean  Survey  (NOS)  charts  were  used  in  the  mesh  devel¬ 
opment:  12222,  Chesapeake  Bay,  Cape  Charles  to  Norfolk  Harbor,  23  January 
1982;  and  12248,  James  River,  Newport  News  to  Jamestown  Island,  4  April  1981. 
These  charts  were  also  used  to  determine  appropriate  depths  for  each  node. 
These  depths  were  not  always  exact  depths  at  nodes  since  some  of  the  elements 
covered  large  variations  in  bathymetry.  In  these  cases,  a  depth  was  chosen 
that  provided  the  proper  cross-sectional  area. 

50.  The  coordinates  of  all  nodes  were  based  on  the  Virginia  State  Grid 
System.  The  grid  is  located  on  the  NOS  charts. 

Development  of  boundary  data 

51.  The  boundary  control  data  for  this  model  were  developed  from  data 
collected  in  the  Chesapeake  Bay  physical  model  for  the  Norfolk  Harbor  and 
channels  deepening  study  (Richards  and  Morton  1983).  There  were  four  test 


conditions  in  the  physical  model: 

Tidal 

Condition  Discharge,  cfs  Range,  ft 

1  200,000  4.8 

2  200,000  3.0 

3  70,000  4.8 

4  70,000  3.0 


52.  The  discharges  shown  in  the  preceding  tabulation  were  the  total 
high  (200,000  cfs)  and  average  (70,000  cfs)  freshwater  discharges  into 
Chesapeake  Bay  from  all  tributaries.  The  tidal  range  was  measured  at  the 
Atlantic  Ocean  physical  model  control  gage  for  spring  (4.8  ft)  and  neap 
(3.0  ft). 

53.  During  the  high  flow,  the  discharge  from  the  James  River  was 
25,363  cfs  and  during  the  average  flow,  it  was  8,877  cfs. 

54.  Figure  16  (tide  stations)  and  Figure  17  (velocity  stations)  show 
the  locations  where  physical  model  data  were  collected  for  the  four  test  con¬ 
ditions.  Each  location  is  described  by  a  range  designator  and  a  station 
number.  The  physical  model  data  were  also  used  to  verify  the  numerical  model 
predictions . 

55.  The  boundary  conditions  needed  were  velocity  data  at  the  upstream 


General  sedimentation  mesh 


Figure  16.  Steady-state  manual  tide  station 
locations 


Figure  17.  Steady-state  velocity  stations 
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boundary  and  in  the  entrance  to  the  Elizabeth  River  and  surface  elevation  data 
for  the  James  River  entrance.  The  velocity  data  had  to  be  depth-integrated 
before  they  could  be  used  as  boundary  data.  These  physical  model  data  pro¬ 
vided  dynamic  updates  hourly  (0-24)  for  the  numerical  model. 

Comparison  of  hydrodynamic  model  results 

56.  Data  from  three  physical  model  stations  were  compared  with  the 
numerical  model  results.  Velocity  data  were  collected  at  all  three  stations 
and  tidal  data  at  one  station.  These  comparisons  were  needed  to  ensure  that 
the  numerical  model  was  calculating  an  accurate  representation  of  the  verti¬ 
cally  averaged  currents  and  tidal  elevations.  Plates  1-4  compare  the  physical 
and  numerical  model  results  for  the  four  conditions  using  data  from  stations 
shown  in  Figures  16  and  17. 

57.  The  comparisons  confirm  the  numerical  model's  ability  to  produce 
proper  currents  and  surface  elevations.  The  comparisons  also  indicate  that 
the  phasing  of  the  tides  during  both  conditions  is  excellent.  The  magnitude 
of  velocity  also  shows  very  good  agreement.  Results  confirm  the  validity  of 
depths  selected  and  the  Manning's  n  values  chosen.  These  comparisons  also 
confirm  the  decision  not  to  include  a  salinity  component  to  this  study 
although  prototype  data  furnished  by  Nichols  (1972)  showed  a  longitudinal 
salinity  gradient.  If  a  salinity  component  had  been  included  in  this  modeling 
effort,  results  would  not  have  been  significantly  improved. 

Hydrodynamic  plan  conditions 

58.  To  simulate  the  proposed  bridge  plan,  only  a  few  changes  in  rough¬ 
ness  were  needed  in  those  areas  that  represent  the  bridge  and  islands.  This 
change  caused  the  two  islands  to  block  the  flow  of  water  and  the  bridge  to 
restrict  the  flow  of  water.  The  boundary  conditions  remained  the  same  for  the 
plan  tests. 

Sediment  transport  modeling 

59.  The  modeling  of  sediment  transport  began  once  the  hydrodynamic  sim¬ 
ulations  were  considered  reasonable.  The  computational  mesh  was  identical 
with  that  of  the  hydrodynamic  model. 

60.  The  following  model  parameters  were  established  during 
verification: 

Parameter 


Time-step  length,  sec 
Crank-Nicholson  implicitness  factor 


Value _ 

1,800  (0.5  hr) 
0.70 


_ Parameter _  Value 

Manning's  n  roughness  0.017 

2 

Critical  shear  stress  for  deposition,  N/m  0.02 

2 

Critical  shear  stress  for  erosion,  N/m 

Top  layer  0.02 

Second  layer  0.06 

2 

Particle  erosion  constant,  kg/m  /sec  0.0012 

Effective  settling  velocity, 

James  River,  m/sec  0.0008 

2 

Dry  weight  of  deposit,  kg/m  500 


The  following  sediment  concentration  boundary  conditions  were  established: 


James  River  (upstream)* 
James  River  (downstream)** 


High  Freshwater 

3 

Inflow,  kg/m 
0.050 
0.10 


Mean  Freshwater 

3 

Inflow,  kg/m 
0.035 
0.010 


*  Values  based  on  Onishi  and  Wise  (1978). 

**  Values  based  on  Nichols  (1972). 

61.  The  sediment  transport  model  parameters  were  adjusted  until  reason¬ 
able  average  shoaling  rates  were  calculated  with  the  concentration  within  a 
reasonable  range.  The  information  in  Nichols  (1972)  was  used  to  determine  an 
average  shoaling  rate  in  the  area  of  interest. 

62.  In  general  Nichols  (1972)  found  the  sedimentation  rates  in  the 
study  area  (excluding  the  navigation  channel)  to  be  less  than  0.07  ft/year. 
After  minor  adjustments,  the  cohesive  (clay)  run  of  the  numerical  sediment 
model,  STUDH,  agreed  well  with  this  observed  shoaling  rate. 

63.  The  sediment  model  was  run  in  the  extrapolation  mode.  This  allowed 
the  model  to  follow  a  sequence  of  a  tidal  cycle  run,  then  extrapolation  and 
another  tidal  cycle.  This  sequence  continued  until  the  bed  had  stabilized  and 
the  model  was  shoaling  at  a  rate  equivalent  to  the  observed  representative 
shoaling  rate. 

64.  It  was  reasonable  to  assume  this  representative  shoaling  value  for 
the  entire  area  of  study  since  Hampton  Roads  is  an  area  of  low  shoaling  and 

no  additional  detailed  data  were  available.  Comparisons  were  made  at  selected 
stations  within  the  actual  area  of  interest.  These  comparisons  confirmed  the 
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ability  of  the  sediment  model  to  adequately  calculate  the  shoaling  in  this 
area  with  adequate  sediment  concentrations. 

65.  Once  the  sediment  model  was  considered  to  be  properly  adjusted,  all 
four  base  and  plan  sets  of  boundary  conditions  were  run.  A  shoaling  index 
(plan  elemental  shoaling  volume  divided  by  base  elemental  shoaling  volume)  was 
used  to  quantify  the  changes  associated  with  the  bridge-tunnel  complex. 
Combining  individual  events 

66.  In  order  to  combine  the  shoaling  indices  for  each  testing  condi¬ 
tion,  as  described  in  paragraph  51,  into  one  yearly  index,  the  duration  of 
each  condition  was  expressed  as  percent  of  a  year: 


Duration 


Condition 

Tide 

Discharge 

percent 

1 

Spring 

High 

1.6 

2 

Neap 

High 

8.4 

3 

Spring 

Mean 

14.4 

4 

Neap 

Mean 

75.6 

Total  100.0 

The  detailed  discussion  of  the  calculations  made  in  establishing  these  per¬ 
centages  is  given  in  Berger  et  al.  (1985). 
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PART  IV:  RESULTS 


Circulation 


Physical  model  data 

67.  The  analysis  of  the  physical  model  data  (Fang  et  al.  1972)  indicated 
there  would  be  no  significant  changes  in  water  levels  or  currents  from  the 
1-664  crossing.  Tidal  elevations  showed  very  little  change.  The  average 
variation  in  the  tidal  range  of  2.5  ft  was  0.1  ft.  There  were  very  subtle 
changes  (less  than  0.5  ppt)  in  the  overall  salinity  structure  in  the  Hampton 
Roads  area  near  Craney  Island.  The  current  regime  was  altered  around  Newport 
News  Point.  The  maximum  velocities  on  both  ebb  and  flood  were  increased  about 
0.75  fps  due  to  the  constriction. 

Surface  current  photographs 

68.  Analysis  of  the  surface  current  photographs  and  the  physical  model 
velocity  data  isolated  those  areas  most  affected  by  the  bridge-tunnel  complex. 
In  general,  the  complex  reduced  the  cross-sectional  area,  causing  velocities 
to  Increase.  In  the  navigation  channel,  changes  in  velocity  were  noted.  The 
north  island  completely  blocked  the  nearshore  flow  and  eddies  developed  both 
upstream  and  downstream  of  the  island.  Strong  currents  developed  on  the  front 
face  of  the  north  island.  Currents  around  the  south  island  Increased  at  the 
north  and  south  ends  of  the  south  island  and  near  the  east  and  west  edges. 

Along  the  bridge,  velocities  increased. 

Navigation  Channel  Sedimentation 

Numerical  model  hydrodynamic  results 

69.  The  hydrodynamic  model,  RMA-2V,  was  run  for  base  (existing)  and 
plan  (proposed  1-664  Bridge-Tunnel  complex  installed)  with  the  four  test  con¬ 
ditions.  Tables  2  and  3  show  the  maximum  ebb  and  flood  depth-averaged  cur¬ 
rents,  respectively.  Locations  of  these  nodes  are  shown  in  Figure  18. 

70.  Boundary  condition  1  had  the  largest  impact  on  navigation  channel 
circulation  of  the  four  boundary  conditions  tested.  Maximum  ebb  velocity  in¬ 
creased  to  about  0.5  f ps  and  maximum  flood  about  0.3  fps  at  node  169,  located 
in  the  channel  between  the  tunnel  islands.  At  node  137,  upstream  of  the  plan, 
the  maximum  ebb  velocity  Increased  about  0.1  fps  and  the  maximum  flood  about 
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0.6  f ps .  At  node  197,  downstream  of  the  plan,  maximum  ebb  velocity  increased 
about  0.6  fps  and  the  maximum  flood  decreased  about  0.1  fps.  Farther  down¬ 
stream  from  the  plan  there  was  no  noticeable  difference  in  the  ebb/flood 
velocity. 

71.  In  general,  downstream  at  node  197,  the  change  in  ebb  velocity  was 
larger  than  that  in  the  flood  velocity.  Upstream  at  node  137,  the  pattern  was 
reversed  with  the  flood  velocity  increasing  more  than  the  ebb  velocity.  In 
the  plan  center  line,  node  169,  the  ebb/flood  changes  were  close  in  magnitude. 
Numerical  model  sedimentation  results 

72.  Thus,  it  can  be  seen  that  velocities  increased  in  the  constriction 
and  in  the  jet  formed  by  the  constriction,  but  the  impact  on  currents  away 
from  the  bridge  was  negligible.  Analysis  of  the  results  indicated  that  the 
Newport  News  Channel  will  remain  basically  stable  with  an  overall  shoaling 
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rate  of  about  0.01  ft/year,  which  represents  a  reduction  from  the  base  condi¬ 
tion  of  0.04  ft/year.  This  shoaling  reduction  was  expected  since  the  hydrody¬ 
namic  numerical  results  indicated  an  average  increase  in  flood  velocities  of 
0.3  fps  and  0.4  fps  in  ebb  velocities. 

73.  The  STUDH  cohesive  sediment  results  from  the  mesh  including  the 
1-664  Bridge-Tunnel  complex  were  very  carefully  compared  with  the  Norfolk 
Harbor  and  channels  deepening  study  plan  runs,  and  no  change  was  noted  in  the 
Norfolk  Harbor  and  Elizabeth  River  Channels. 

General  Sedimentation 


Numerical  model  hydrodynamic  results 

74.  Plates  5-12  are  vector  plots  of  the  vertically  averaged  velocities 
computed  at  the  computation  locations  in  the  area  of  interest  (entire  computa¬ 
tional  mesh  shown  in  Figure  15).  These  plots  show  maximum  ebb  (hour  15)  and 
flood  (hour  21)  velocities  for  the  four  base  and  plan  conditions  that  repre¬ 
sent  various  combinations  of  tide  ranges  (spring  and  neap)  and  discharge  con¬ 
ditions  (high  and  mean) . 

75.  The  bridge  caused  maximum  velocities  to  increase  downstream  on  ebb 
and  upstream  on  flood.  These  increases  were  small  and  would  have  a  minor 
impact  on  circulation.  Between  the  two  islands,  both  upstream  and  downstream 
velocities  were  increased.  This  effect  diminished  as  the  flow  moved  toward 
Craney  Island.  Directly  in  front  of  Craney  Island,  the  flow  reestablished 
itself  and  was  essentially  unchanged  from  the  base  test. 

Numerical  model  sedimentation  results 

76.  Figure  19  shows  the  yearly  simulated  computed  shoaling  indices  for 
selected  elements.  These  elements  were  selected  to  evaluate  sedimentation 
changes  within  the  confines  of  the  oyster  beds. 

77.  The  sedimentation  on  oyster  bed  locations  6  and  7  (Figure  14), 
located  on  either  side  of  the  proposed  bridge,  remained  unchanged  or  decreased 
slightly.  This  was  as  expected  since  velocities  were  slightly  increased  in 
these  areas. 

78.  Sedimentation  remained  unchanged  on  oyster  beds  2  and  3  (Fig¬ 
ure  14).  This  area  is  effectively  out  of  the  influence  of  the  1-664  Bridge- 
Tunnel  complex. 

79.  Sedimentation  did  not  change  in  the  majority  of  the  public  clamming 
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grounds  (shown  as  Area  1  in  Figure  14).  However,  sedimentation  increased  near 
the  north  island  as  a  result  of  nearshore  circulation  blocked  by  the  island 
and  generation  of  an  eddy.  This  eddy  could  be  moving  sediment  into  this  area. 
The  prototype  sedimentation  in  this  area  is  relatively  minor,  and  the  in¬ 
creases  caused  by  the  north  island  would  probably  be  local. 
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PART  V:  CONCLUSIONS 


80.  The  TABS-2  model  of  the  lower  James  River  accurately  reproduced 
currents  and  surface  elevations.  The  phasing  of  the  tides  was  excellent  and 
the  magnitude  of  velocity  also  showed  very  good  agreement.  It  can  be  expected 
to  provide  reliable  guidance  for  hydrodynamic  and  sedimentation  impacts  of 
proposed  construction. 

81.  The  model  test  results  showed  that  the  1-664  Bridge-Tunnel  had  the 
following  effects: 

a.  Circulation  changes  and  resulting  flushing  changes  were 
localized  with  minimal  effects  on  the  general  circulation  of 
the  lower  James  River. 

b.  Sedimentation  was  generally  unchanged  or  reduced  except  on 
either  side  of  the  north  island  where  increases  will  be 
expected.  The  areas  experiencing  unchanged  or  slightly  reduced 
sedimentation  rates  included  the  oyster  grounds,  the  Elizabeth 
River  and  Norfolk  Harbor  Channels,  and  the  Newport  News 
Channel . 
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Table  1 

Velocity  and  Salinity  Data  Collected 


During  Circulation  Study  by  VIMS 


Station 


Depth 


_ Base _ 

Velocity  Salini 


Veloci 


*  Indicates  data  collected  at  this  location 


r 

r 


[« 

1 


n 


fi 


Table  1  (Concluded) 


Base 

Plan 

1A 

Station 

Depth 

Velocity 

Salinity 

Velocity 

Salinity 

4A 

S 

M 

* 

* 

* 

* 

B 

5A 

S 

* 

* 

M 

* 

* 

B 

* 

* 

5B 

S 

* 

* 

* 

* 

M 

* 

* 

B 

* 

* 

* 

* 

5C 

S 

* 

* 

M 

* 

* 

B 

* 

* 

5D 

S 

M 

* 

* 

* 

* 

B 

5H 

S 

* 

* 

* 

* 

M 

* 

* 

B 

* 

* 

* 

* 

6A 

S 

* 

* 

* 

* 

M 

* 

* 

B 

* 

* 

* 

* 

6B 

S 

M 

* 

* 

* 

* 

B 

6C 

S 

M 

* 

* 

* 

* 

B 

6D 

S 

M 

* 

* 

B 

7 

S 

* 

* 

M 

* 

* 

* 

* 

B 

* 

* 

*  Indicates  data  collected  at  this  location. 


Table  2 

Maximum  Ebb  Velocities 
Test  Conditions 


High  Discharge,  fps _  _ Mean  Discharge ,  fps 


Spring 

Tide 

Neap  Tide 

Spring  Tide 

Neap  Tide 

Node 

Base  1 

Plan  1 

Base  2 

Plan  2 

Base  3 

Plan  3 

Base  4 

Plan  4 

137 

1.94 

1.04 

2.17 

2.28 

1.57 

1.63 

169 

2.27 

2.79 

1.12 

1.51 

2.50 

3.03 

1.67 

2.11 

197 

2.68 

3.32 

1.37 

1.77 

2.88 

3.57 

2.00 

2.49 

223 

1.80 

1.82 

0.98 

0.99 

2.01 

2.09 

1.39 

1.27 

262 

1.75 

1.70 

0.93 

1.88 

1.80 

1.22 

1.20 

Table  3 

Maximum  Flood  Velocities 
Test  Conditions 


High  Discharge,  fps 

Mean  Discharge ,  fps 

Spring 

Tide 

Neap 

Tide 

Spring  Tide 

Neap  Tide 

Node 

Base  1 

Plan  1 

Base  2 

Plan  2 

Base  3 

Plan  3 

Base  4 

Plan  4 

137 

1.51 

2.09 

1.03 

1.50 

1.47 

2.17 

1.23 

1.70 

169 

2.40 

2.74 

1.61 

1.90 

2.27 

2.65 

1.85 

2.16 

197 

1.90 

1.78 

1.30 

1.28 

1.75 

1.73 

1.46 

1.39 

223 

1.50 

1.44 

1.20 

1.18 

1.59 

1.57 

1.13 

1.09 

262 

0.97 

0.92 

0.85 

0.81 

1.19 

1.16 

0.73 

0.70 

NOTE:  Base  1  Indicates  the  base  configuration  with  boundary  condition  1; 
Plan  1  indicates  the  bridge/tunnel  configuration  with  boundary 
condition  1;  etc. 
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